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By R. Betchov 


We study here the properties of a hot-wire anemometer under the 
supposition that the heat transfer from the wire to the air depends, 
first, on the difference in temperature and, second, on the square of 
that difference. This latter hypothesis is confirmed by experience, 
and the consequences might be of great importance, that effect of non- 
linearity is stronger than the effect of thermal conduction. 


I. THE NONLINEAR LAW OF KING 

The heat quantity Q removed per second by an air stream V from 
a wire of the diameter d and unit length is given by King in the form 

Q = (a + by/Vd)T (l) 

where T denotes the temperature difference between wire and air. King 1 
calculation, approximately confirmed by experience, yields 

a = % x b = \Z2rtX f S f c T (2) 


with x 1 = thermal conductivity of the air, 5* = density of the air, 
and c T = specific heat of the air for constant volume. 

These quantities may vary with T, and experience shows that a 
increases while b remains practically constant. Intuitively, one may 
interpret this effect by saying that the air in contact with the wire 
is heated which increases its conductivity. In compensation, its density 
decreases because the pressure varies only very slightly. Obviously, 
the effects on x* and 5* compensate one another, and only a varies. 


*"Theorie non-lineaire de l’anemometre a fil chaud. n Koninklijke 
Nederlandsche Akademie van Wetenschappen. Mededeling No. 6l uit het 
Laboratorium voor Aero- en Hydrodynamica der Technische Hoge school te 
Delft. Reprinted from Proceedings Vol. LII, No. 3, 19^9, pp. 195-207. 



2 


NACA TM 1346 


King states in his original report (ref. 1 ) that a increases by 
0.114 percent per degree; he also describes there an effect of the 
diameter on that term a which we shall not discuss here. Thus it is 
advisable to write 



(3) 


where 7 takes the nonlinearity into account. 

One must not forget the hypotheses on which King bases his calcu- 
lation: he contends that the air flow is without viscosity and that 

the heat flow in the immediate proximity of the wire is constant. He 
uses the specific heat at constant volume although the pressure is 
certainly more constant than the density. For that reason, we consider 
equation (3) as an empirical relation, valid for the wire unit length, 
and would wish to see King ! s problem made the subject of a more thorough 
investigation. 

Here we intend to study the effect of the term 7 on the properties 
of the hot wire; we simplify the notation by introducing P so that 

P = |\/vd (4) 


Q = a(l + P + 7T)T ( 5 ) 


II. GENERAL EQUATION OF THE HOT WIRE 

We shall use the following symbols: 

S resistance of the wire, per unit length, at operating 

temperature 

S Q resistance of the wire, per unit length, at ambient temperature 

I intensity of the electric current heating the wire 

a coefficient of the variation of S according to the 

temperature 

m weight of the wire, per unit length 
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c specific heat of the wire, in joule/gram and degree 

% coefficient of the thermal conductivity of the wire in watt/cm 

and degree 

a wire cross-sectional area 

Z semilength of the wire 

t time 

x coordinate of position, varying from Z to -Z 

We put 


A = 


ccS, 


-(1 + P) = 


a + 


b\/vd 


aS r 


(6) 


The equation of the hot wire must express the equilibrium between 
the heat supplied per second, the heat removed by the air stream, the 
heat required to raise the temperature of the wire, and the heat trans- 
mitted by conduction. One obtains 


SI 2 = A(S - S 0 ) + -^-(S - S 0 ) 2 + S— 
«2s - 2 aSr 


dS 


v rr 


*2 C 


For the steady-state case, and introducing the parameters 




Z = 



(S - S 0 ) 




2 ay I 2 /A 2 

3 a 2s 0 (1 - I 2 /A)2 


2 7 1 I 2 /A 

3 a 1 + P (l _ i2/a)2 


(8) 


one obtains the equation (7) in the form 


Z + | GZ 2 



(9) 
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III. EXACT INTEGRATION OF THE STATIC CASE 

Multiplying equation (9) by dz/dy and integrating, one obtains, 
with a constant 


— = \/gZ 3 + Z 2 - 2Z + const (10) 

dy 


One notes that dz/c)y is zero for a negative value of Z and may 
be zero for two positive values of Z. At the ends of the wire, one 
has S = S Q and Z = 0; at the center, Z must be positive and dz/dy 
zero. The range of interest for us lies, therefore, between Z = 0 and 

the first positive root which gives — = 0 which we shall denote by 

oy 


Z = B. We put 


z(y) = B - X 2 (y) 


( 11 ) 


By virtue of the relation 


qb3 + b 2 - 2B + const = 0 (12) 

one obtains 

4(dx/dy) 2 = -GX^ + EX 2 + D (13) 


with 


E = 1 + 3GB D = 2(1 - B) - 3GB 2 (l4) 

Following, we shall consider B as a new integration constant 
indicating the temperature in the middle of the wire. At the center of 
the wire, one has X = 0 and 


4(dX/dy) 2 = D 


(15) 
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which shows that D is positive and generally small. At the ends of 
the wire, Z = 0 and X = ±/b. 

The roots of equation (13) are 


x 2 = -E * VE 2 + 4GD 
-2G 

Introducing the parameter (3 so that 


(16) 


sinh (3 = sh p = — 
E 


(17) 


one can write equation (13) in the form 


4(8x/6y) 2 = G ~(ch 3 


+ 1 ) - x c 


f(ch P - 1) + X : 


(18) 


We define the angle cp, function of y, so that 


=(-> 


sin cp 


E ch (3 ./ o p - 
]/l - k^sin^cp 


(19) 


with 


k 2 = 


ch p + 1 
2 ch p 


Equation (18) then he comes 


dtp 


\fl - k 2 sin^p 


= We ch |3 dy 


( 20 ) 


( 21 ) 
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and we obtain the elliptic integral of the first kind 


U(k>cp) = = | }/i ch Py 

-'o \Ji ■ ° 2 


k 2 sin^cp 


( 22 ) 


The variable y varies from -I to I , with 


I = l/V 


(23) 


Q 

and we have, for y = I and = B 


» r$. max 

-VEchl = / 

2 Jn 


dcp 


0 \f± - k 2 sin 2 cp 


(24) 


with 


Vb 


'Vi 


sin cp 


max 


E ch P \Jl - k 2 sin^p 


max 


(25) 


This last equation may be written 


1 - v2 


a k^ + 


smV EB ch P 


nnax 


(26) 


From equation (20) one may deduce 


— = 1 + tanh2(p/2) 


( 27 ) 


With 3 ranging from 0 to +°°, k 2 varies from 1 to 0.5. 
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With G and B known, one may calculate successively D, E, (3, 
k, and <p max . A table of U(k,cp) then permits us to calculate 

Figure 1 gives the results obtained by this procedure and allows - 
starting out from a prescribed wire and with 7 known - to deter min e 
B from G and I . 

The temperature distribution over the wire is given by Z as a 
function of y, or by X as a function of cp and cp a function of y. 
The relation cp(y) is given by the quotient of equations (22) and (24), 
namely 


y 

£ 


r 9 &p 

w p \/l - k^sin^cp 


L 


^max 

0 \J± - k2 S in2cp 


dcp 


(28) 


From equations (19) and (25) , one obtains a relation between cp and X, 
namely 


x£ _ sin^cp 1 - k sin^cp max 
B sin2c Pmax 1 “ k2sin2< ? 


(29) 


The total resistance R of the wire is given by 



S dx = 


2S 0 I 2 
A - I 2 



dy + 2S 0 Z 


(30) 


We introduce the cold resistance 


R 0 and the function X 


R 




X2 dy 


0 


(31) 
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We replace X according to equation (19) and dy according to equa- 
tion (21), namely 


r r 

o 


-,9- 


max 


* - * ‘ fTP'Z - J 0 ~ <32) 


The integral is equal to 2<3u/<3k 2 arid we obtain 


R 


_ R o = 3^L/ B _ 

° A - I2\ E ch 3 U / 


(33) 


The values of 8u/<3k 2 can be deduced from a good table of U(k,cp) 
with sufficient approximation. 

If the wire were perfect, the expression in parentheses in equa- 
tion (33) would have to be replaced by unity} therefore, we shall intro- 
duce the quantity M so that 


M = 1 - B + gI) au / Skg 
E ch p U 


The formula (33) then gives us 


R = 


P 1 - mi 2 /a 
E ° 


(34) 


(35) 


and the important relation 


RI^ 

R - R, 


= A 1 + 


M 

1 - M 


(1 



(36) 


This last equation permits easy determination of the wire charac- 
teristics because R, R Q , and I can be measured accurately and because 

the curves obtained as functions of R/R q , for instance, indicate directly 
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the effects of conduction and of nonlinearity. We calculated the values 

of — - — as a function of G and £1 figure 2 shows our results. 

1 - M 

A good approximation is given by 


1 




(37) 


B 0 


+ y/l + 6 g - l 

3G 


(38) 


IV. A FEW USEFUL APPROXIMATIONS 


In performing the calculations necessary for the plotting of fig- 
ure 1, we have noted that one may assign to k the value unity without 
introducing large errors. 

This implies (3 = 0 and equation (19) then gives 


X 



tan cp 


(39) 


The integral (22) becomes 


U = 




dcp 


Jo cos <P 


1 L 1 ± sin cp = 1 iEy 

2 1 - sin cp 2 


(40) 


whence, one deduces 


ch(2U) - 1 * 3ln ^ = ch(VEy) 

1 - sin^cp 


( 41 ) 
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X 2 = ^(ch \f¥y - 1) 


(42) 


At the limits, one has 


B = — (ch (Ei - 1) 
2E 


(43) 


which gives 


B L _ ch \[Ey \ 

1 - l/ch \[Et \ ch {Eij 


(44) 


In order to calculate M, one must put 


lim Su/Sk 2 = i r ^ dtp = i/'iiSJE 

k — >1 


2 J 0 cos ’ 4 co s 2 cp 


- U(cp) 


(45) 


which gives 


M = 1 - 


B 


Th \/||\ 

1 - l/ch /Ell \/E| ) 


1 - 


(46) 


One can see that, due to the nonlinearity, 1 is replaced by v^El, 
and the central temperature is lowered. 

If one takes equation (44) as solution of equation ( 9 ), one sees 
that the equation is satisfied for a term of approximately (ch \/Ey/ch \/e |) 2 , 
and that B is given approximately by 


B = ^V 1 + ^j G - ^ (1 _ l/ch m) (47) 
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with 


E - \/l + 6 G 


(43) 


We shall take an example that represents an extreme case. We 
choose a platinum wire with 10 percent of iridium, a diameter of 7 microns, 
and a length of 22 = 1.14 mm. Exposed to an air stream of 5 m/sec and 
heated with 75 mA, it gives us 


P = 2.9 A = 1.2 X 10 -2 I 2 /A = 0.47 

2 * = 0.15 mm t = 3-75 

Assuming y = 1.2 X 10“3 and with the aid of figures 1 and 2, we 
determined 


G = 0.25 B = 0.76 E = I .56 M = 0.4 


The other parameters have the following calculated values: 


k = 0.99756 


^max 


=79.5° 


sh (3 = 0.14 


D = 0.0477 U = 2.34 


In figure 3, we show the profile of the temperatures calculated 
exactly, the profile calculated with the approximation k = 1 , and the 
profile calculated with 7=0. It can be seen that the nonlinearity 
offers a more uniform temperature distribution, and that the approxi- 
mation is sufficient. 

By means of equation (35) one calculates 5 — = 1.5, whereas the 

n 0 

calculation with 7=0 would give the result 1.9. The mean tempera- 

ture giving — = 1.5 would be 38 O 0 . 

R o 

As to the term RI 2 /R - R 0 , it changes from the value 1.24 A when 
the current is very weak to the value 1.35 A when I attains 75 mA. 
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R R 

It varies therefore by about 10 percent between — = 1 and — = 1.5 

Ro Ro 

which is of the same order as the variations observed. 

Thus this magnitude is constant at 10 percent and King's law is 
verified; however, we shall see later on that the thermal inertia is 
very different from the expected value. 


V. THE DYNAMIC EQUATION OF THE HOT WIRE 


In order to calculate the variation of the total resistance of the 
wire when the current or the air stream fluctuate, one must go back to 
the equation (7). Replacing in this equation I, S, and V (contained 

in the term A) by I + ie'Prt, s + , and V + ve^ 0 ^, one obtains 

after suppressing the terms of the order zero, two, and more as well as 
the factor eJ^ 


2SI1 + I 2 s = -SL. I( S - s 0 ) + As + 
2 aS Q V 


2ay (S - S Q )s + 

a 2 s/ 


me 



jeos 


xa <3^s 


(49) 


We introduce 



(50) 


and, identical to the formula (9) of our Mededeling'No. 55 (ref. 5) 


cn* = — —(A - I 2 ) 
me 

One then obtains, after introduction of Z 


(51) 


o i , o i I 2 /A „ 1 P Z v 

1 1 1 - I 2 / A 21+ ?1- i2/ A V 


z(l + 3GZ + ) 


6y 2 


( 52 ) 
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The function Z(y) appears twice in this differential equation, and we 
must utilize the approximation k = 1 in order to avoid great compli- 
cations. Writing the expression Z according to equation (44), one 
obtains 


a 2 z L 

+ z 1 + 

6y 2 V 


3GB 

1 - l/ch VEl 


jcjo/cu* 


3GB 

ch \/E| - 1 


ch v/Ey 


“ Ci 


ch y/Ey 
2 ch v/E| 


(53) 


with 


Cl . 2 i/l + 4_J| . i Ip i 

I\ 1 - I 2 /A 1 - 1 /ch /E$/ 2 V\1 + P x _ j2/ a 1 - l/ch \/E£/ 

(54) 

c = _2 i i2 / a B + 1 v / P 1 B ^ 

2 1 1 - I 2 /A 1 - l/ch v/E| 2 V^l + P -l _ j2/ a 1 - l/ch \/Eiy 

(55) 


The solution without the second member of equation (53) is a 
Mathieu function taken for a purely imaginary value of the argument-, 
but if \/e| is large enough, for instance , more than 4, one may neglect 
the term with ch \/Ey of the member on the left side of equation (53)* 
Actually it is important only when y is close to however, we shall 
take as a limiting condition z(|) = 0, and the product z ch VEy will 
never be important. We shall also neglect the term l/ch \/e£ in the 
denominator of one of the terms on the left side of equation (53) which 
amounts to taking Z = B in the factor term of z. 

• Taking the definition of E into account, one then has 


"~§ + z( E + W®*) = Cl + c 2 ( 56 ) 

dy 2 ch VE£ 

the solution of which - null for y = ±1 - is 

C 1 A _ ch p/Ey \ + C 2 /ch \[Ey _ ch_p\/Ey\ 

Ep2 \ ch V m) jo)/o)*\ch \[Ei ch p\fEt) 


z 


(57) 
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with 


p = \/l + jjo/Eco* 


(58) 


It can be seen that due to the nonlinearity the characteristic 
frequency u# is replaced by a new frequency which is higher by a fac- 
tor E. (Compare with equation (8) of ref. 5.) The amplitudes according 
to Ci and C2 also are modified by the nonlinearity. 

The variation r of the total resistance R will be 



(59) 


The integration gives 



A - I 2 


tanh p _^2 W tanh ^e| . tanh pi/E| 


r>\kt 


We assumed above that v^El is sufficiently large; also, we may 
assign to the hyperbolic tangents the value unity (the presence of a 
complex argument is here not of importance). 

One then obtains 


r 



%>(1 - 1 / p ^ E !) + 

Ep^ 


C 2 1 ( V -1 \1 

jo)/o# y^EiV p /[ 


(61) 


VI. RESPONSE TO A FLUCTUATION OF THE CURRENT 

If one assumes v = 0 in equations (54) and (55), one may trans- 
form equation (6l), suppressing the terms containing l/ch /e? 


2iR Q I 2 /A 

(1 - I 2 /A(l - B) ) (l - l/p/li) bi 2 /a O)*; 


(1 - i 2 /a) 2 

Ep2 - m >1 

( p t 


rl = 


(62) 



NACA TM 1346 


15 


This formula gives the alternating electromotive force produced 
at the boundaries of the hot wire by the modulation current i in 
addition to the normal electromotive force Ri. Although this formula 
seems to be complicated, it can be adapted to the needs of practice. 

If the frequency cd/2k tends toward infinity, equation ( 62 ) becomes 


= 2 iR o l2 / A (1 - I 2 /a(1 - B) - I 2 /AB/y/Eg) 

(1 - i 2 /a ) 2 > 


If one takes into account that equation (35) gives R, that equation (46) 
gives M, and that equation (51) gives u)*, one finds 


ri = -^r- iCRI 2 
(x>/2n 


with 



jtmc 


(64) ( 65 ) 


The constant C corresponds to that of our former publications, 
and equation (64) shows that the electromotive force ri taken at high 
frequency permits to measure C without being impeded by either con- 
duction or nonlinearity. The method described in reference 4 is there- 
fore indicated rather than the one consisting of measuring the phase 
displacements, with od being of the order of ay*. 


One may immediately verify this point by assuming in the equa- 
tion (49) as very large, thus making the effect of the terms of conduction 
and of nonlinearity negligible. 


When a) tends toward zero, the electromotive force becomes 


ri = 


21RqI 2 /A j f 

(1 - I2/A)2 E 1 




( 66 ) 


Thus the complex function ri according to equation (62) will change 
from equation (64) into equation (66) when oj varies from 0 to a large 
value. The complex trace of this function gives practically a semicircle 
which permits to put approximately 


rI = rl(o> = 0) 

1 + jco/o )** 


( 67 ) 
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where a#* denotes the effective characteristic frequency. In fig- 
ure 4 , we plotted the semicircle and indicated a few values of tu/cu**. 
From equation ( 62 ), and in the case of the preceding example, we calcu- 
lated the electromotive forces rl for a few values of cu/Ecu*. One 
can see that the two functions blend, at low frequencies, if one assumes 
to** = l.lEtu* = 1.7 &d*. At high frequency, equations ( 64 ) and (67) will 
be equal which permits calculation of a satisfactory value of cu** when 
the effect of thermal inertia is important. One obtains 


CD** 


Eco* 


1 - l/v/E£ 


1 - r 


- j 


1 - i 2 /am 


B 


( 68 ) 


In the case treated one finds oj** = 1.23&U*, that is, cu** = 1.93^*: 
the effective characteristic frequency is almost twice the expected 
value) therefore, the approximation (67) gives a correct plot of the 
function, but the phases according to equation ( 68 ) will be only within 
a 10 -percent accuracy. 

The denominator of equation ( 68 ) depends chiefly on I, and it 
increases the characteristic frequency. Instead of compensating each 
other as in the static case, the two effects reduce the thermal inertia. 

Intuitively, one may say that the conduction shortens the hot part 
of the wire and thus reduces the heat required for modifying the central 
temperaturej the nonlinearity depends on the presence of hot air around 
the wire, and the thermal inertia of the air is negligible which improves 
the spherical response of the anemometer. 

When the wire is dusty, the quantity of immobile air is greater, 
and the experience shows that the term a in equation ( 3 ) is increased 
while b remains unchanged. One must therefore expect a dynamic action 
of the dust of the wire to the extent that E is modified. The dynamic 
effect may be more important than the static effect. 

The wire in the quoted example demonstrates that, with a term 
RI 2 /R - R 0 constant at 10 percent, the characteristic frequency may be 

almost twice the normally foreseen value. 
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VII. RESPONSE TO A FLUCTUATION OF THE AIR STREAM 


Assuming i = 0 in the formulas (5*0 and (55) , and maintaining v , 
one may transform equation (6l) into 


rl = - — — 


1 v P 


PqI 3 /A 




B 


2 V 1 + P (i _ i2/ A )2 i _ l/ch l/i& 


{ 


1 - l/p v/E| co* p - 1 1 | 

jcu p Vei J 


Ep £ 


(69) 


With cd tending toward zero, one has 


r 


T 1 V P 
rl 


r q i 3 /a 


B 


2 V ! + P (i _ p2/ A )2 p _ p/ ch 


1 - * 


3 _1_ 
2 i/es 


E 


J 


(70) 


If cd tends toward infinity, one has 


T 1 v P 
rl = 


r q i 3 /a 


2 V 1 + P (i _ p2/ A )2 p _ p/ ch y/E| 


1 - 


\/eI 


J(d/cD* 


(71) 


and one may approximately replace equation (69) by the semicircle 


= rl(cp = 0) 

1 + jcd/cd** 


(72) 


with 


CD** = Ecd* 1 

1 - | 1/v/El 


(73) 


which, in the case of the example treated above, gives cd** = I.8I0*. 
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Thus there is, on principle, no equality between the dynamic reac- 
tion to a variation i and to a variation v. 

This difficulty arises due to the term Ve|, namely to the conduction 
the nonlinearity tends toward diminishing its importance (factor \/E). 

We hope to publish some empirical results, and the calculation of 
the differences indicated by different authors, in the near future. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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Figure 2.- Values of M/1 -M according to G and 
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Figure 4.- Study of the tension rl. The semicircle corresponds to 
formula (67); the points give rl, according to (62), for several values 
of u/Ew*. 
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